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Abstract 

In this paper we describe the implementation of the SM process 77 — > 77 through a fermion and 
boson loops into the framework of SANC system. The computations of this process takes into account 
non-zero mass of loop particles. We briefly describe additional precomputation modules used for 
calculation of massive box diagrams. We present the covariant and helicity amplitudes for this 
process, some particular cases of Do and Co Passarino-Veltman functions and also numerical results 
of corresponding SANC module evaluation. Whenever possible, we compare the results with those 
existing in the literature. 
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1 Introduction 



SANC is a computer system for semi-automatic calculations of realistic and pseudo-observables for 
various processes of elementary particle interactions "from the SM Lagrangian to event distributions" at 
the one- loop precision level for the present and future colliders — TEVATRON, LHC, electron Linear 
Colliders (ILC, CLIC), muon factories and others. To learn more about available processes in SANC see 
the description in [IJ [2] and look at our home pages at JINR and CERN [3] . 

Light-by-light scattering is one of the most fundamental processes. It proceeds via one-loop box 
diagrams containing charged particles. The first results for the QED low energy limit of this process 
were obtained by Euler [?]. Then Karplus and Neumann [5] found a solution for QED in general but 
complicated way. The QED cross sections in the high energy limit were calculated by Ahiezer [5]. 
Nowadays there are computations for 77 — > 77 process in the electorweak Standard Model [7l HJ O [10] 
and even for two-loop QCD and QED corrections [IT] . 

In this paper we describe the implementation of the SM process 77 — > 77 through fermion and boson 
loops and corresponding precomputation modules into the framework of SANC system. The computations 
of this process take into account non-zero mass of the loop particles. 

One should emphasize that the obtained building blocks and procedures of precomputation for box 
diagrams in QED and EW (as in 77 — > 77) is the first step in the creation of environment for implemen- 
tation of the similar four-bosons processes in the Standard Model (like 77 — > ZH , 77 — > ZZ [TO]) and in 
QCD (like gg — ► 77, gg — ► ZZ, gg -> W + W~ etc.)- 

The paper is organized as follows: 



First we discuss some notations and common expression for cross section in section 2.1 

In section [2~2| we discuss diagrams for 77 — > 77 process and covariant amplitude tensor structure. 

The helicities amplitudes approach [TJ [T^] and their expressions for light-by- light scattering in general 



(massive) and in limiting (masslcss) cases arc listed in section 2.3 

In section [3] we shortly describe precomputation strategy and the place of this process on the SANC 
process tree. 



The implementation of analitycal results and the SANC modules concept we describe in section 3.2 
At last in section |3.3| one can find the numerical result and comparisons with those existing in the 
literature. 

Additionally, in Appendix section we list answers for particular cases of Dq, Cq and Bq Passarino- 
Veltman (PV) functions [T3] (see also the book [3]), which are needed for calculation of light-by- light 
scattering through massive and massless loop particles. Finally, we present strings and basis for covariant 
amplitude. 

2 Light-by- light scattering process 

2.1 Notation, cross section 

The 4- momenta of incoming photons are denoted by p\ and P2, of the outgoing ones — by p 3 and p^. 
The amplitudes are calculated for scattering of real photons, that is p\ — , p\ — , p\ — , p\ = 0. 
The 4- momentum conservation law reads: 

Pi +P2 ~P3 - Pi = 0. (1) 

The Mandelstam variables areQ 

s = -(P1+P2) 2 = -2pi -p 2 , t= -(pi -p 3 ) 2 = 2px -p 3 , 

u = —(pi — pi) 2 = 2pi ■ pi, s + 1 + u = . (2) 



-'^Note, that in SANC wc use Pauli metric. 
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For the 2 — > 2 77 — > 77 process, the cross section has the form: 



(3) 



where j = 4-y/ (P1P2) 2 is the flux, yl 77 _^ 77 is the covariant amplitude (CA) of the process, and <i<I>( 2 ) is 
the two body phase space: 

dW (p§) d 4 p 4 5 (pi) 



d<S>^ = (2^6 ( Pl + P2 - P3 - P4 ) 
For the differential cross section one gets: 



(2tt) 3 (2tt) 2 
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(5) 



where s = 4oj 2 , u> is the photons energy and 9 — the scattering angle in the center of mass system (cms) . 

2.2 Covariant amplitude 

The covariant one-loop amplitude corresponds to a result of the straightforward standard calculation 
of all diagrams contributing to a given process at Born (tree) and one-loop (1-loop) levels. 

The CA is being represented in a certain basis, made of strings of Dirac matrices and/or external 
momenta (structures), contracted with polarization vectors of vector bosons, e(p), if any. The amplitude 
also contains kinematical factors and coupling constants and is parameterized by a certain number of 
Form Factors (FFs), which are denoted by in general with an index i labeling the corresponding 
structure. The number of T\ is equal to the number of independent structures. 

The 77 — > 77 process in QFT appears due to non-linear effects of interaction with vacuum, so this 
process has no Born or tree level. Corresponding diagrams start from the one-loop level and in QED there 
are box diagrams with four internal fermions of equal mass. The number of not identical diagrams (or 
topologies) is equal to six. But three of them differ from another only by the orientation of the internal 
fermionic loop, giving the same contribution or a factor 2 to the amplitude. So, only three topologies (st, 
su and ut channels) remain which are related by simple permutations of external photons in the diagrams 
shown in Figure [lj su-channcl is obtained from st-channel by p 3 <-> p4 rotation and ut-channel — by 
P2 <-> P3- The sum of these fermionic diagrams is a gauge invariant in each generation of particles. 



A(-p4) 



A(-p3) 



A(-p4) 



A(-p3) 




A(+pl) A(-p4) 




Af+pl) A( + p2) 




A(+pl) 



A(+p2) 

(a) st-channel 



A(+p2) 

(b) su-channcl 



A(-p3) 

(c) ut-channcl 



Figure 1: 77 — > 77 process QED diagrams 

In the EW boson sector we have three types of diagrams to classify: box topologies, pinch topologies 
and fish topologies (shown in Figure [2]). There are three channels of each topologies (st-, su- and ut- 
channel as in QED) and we have W + , W~ , cj) + , $T and X + , X~ (bosons and ghosts) as internal particles 
in i?£-gauge theory. 
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As in fermionic part we can choose only positive or negative charge bosons and X + , X~ ghosts to 
appear as loop particles and multiply the result by factor 2 to dissmiss the double counting diagrams, 
which differ from another only by the orientation of the loop charge flow. 

So, we have three structures (3 channels) in box type of diagrams, twelve structures (3 channels by 4 
corresponding pinches) in pinch type and six structures (3 channels by 2 corresponding combinations of 
propagators — direct and crossed) in fish type of diagrams — each of ones is a sum of the appropriate 
sets of loop particles diagrams. 




The full CA of given process for off-shell photons (ptCi ^ 0) with corresponding combinatorical factors 
can be written as sum of bosonic part minus fermionic and ghost part: 




In terms of Lorenz-structures we have: 

43 



^ 77 ^ 77 = | ^ osons (a , t , u) + ^™ ( s , t , u) T^ v . (7) 

i=\ 



The Ti are normilized by corresponding factors for fermion and boson parts: 

C formions = 8a 2 QjN c , 

^bosons = 12q , 2j (g) 

where a is the fine structure constant, Qf is the fraction of charge of loop fermion in units of electron 
charge e, iV c is the number of colours for given fermion, T° v are tensors defined with an aid of auxiliary 
strings Tj presented in the Appendix section IV The off- shell process contains 43 basis elements, but for 
the on-shell real photons we need only first 10 structures. 
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Thus, one gets a minimal number of tensor structures of the CA. It can be written in an explicit form 
with an aid of scalar Ti. All masses and other parameters dependences are included into these T^ It 
is important to emphasize that each basis elements T^^ LV , i = 1 43 of tensor structure with Lorenz 
indices is four times transversal with respect to each external photon: 

Tf» v p a = Tf» v Vf3 = Tf» v Pll = Tf» v Vv = . (9) 

Ti are scalar coefficients in front of basis structures of the CA — projections of CA to complete basis 
expressions T° " . They are presented as some combinations of scalar Passarino-Veltman functions Aq, 
B , Co, D |13j and depended on invariants s ,t ,u and also on fermions and bosons masses. They do not 
contain UV poles. 

The number of terms in Ti equals to thousands in the case of non-zero mass of the loop particles, but 
this number reduces greatly for zero loop fermion mass. Full answer for Ti one can find in the computer 
system SANC. The client part is available on servers [3]. 



2.3 Helicity amplitudes 

In SANC we use helicity amplitudes approach [TJ [TJ] . 

In the expression for CA as one can see in subsection 2.2 one has tensor structures and a set of scalar 
Ti. To calculate an observable quantity, such as cross section, one needs to make amplitude square, 
calculate, in general, products of Dirac spinors and contract Lorenz indices with polarization vector. In 
the standard approach making amplitude square one gets squares for each diagram and their interferences. 
This leads to a huge number of terms. 

In the helicity amplitudes approach we also derive tensor structure and Ti. But the next step is a 
projection to helicity basis and as a result one gets a set of non-interfering amplitudes, since all of them 
are characterized by different set of helicity quantum numbers. In this way we can distinguish calculations 
of Dirac spinors, if they are needed, and contraction of Lorenz indices from calculations of Ti. We can 
do this before making square of amplitude. So we get a profit on calculation time (less amount of terms 
due to zero interference) and also more clear step- by-step control. 

For the process 77 — > 77 one gets: 
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Note, the total number of HAs for this process is equal to 16. This corresponds to different com- 
binations of external particles spin projections. For 77 — > 77 processes there are 4 photons with two 
polarizations '+' and '— ', so the total number is 2 • 2 ■ 2 ■ 2 = 16. Helicity amplitudes are scalar expressions 
with factors equal to C fcimlons for fermions and C bosons for bosons 

For bosons part we have: 
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The analytical answers for scalar Passarino-Veltman functions A , B , Cq, D |13j are presented in 
corresponding Appendix sections. 

For fermions part we have similar to bosons result with corresponding loop particles masses, but with 

opposite sign and there is a difference in helicities conserving part of amplitude (+ + ++ and H 1 — 

permutations), the non-conserving helicities (+ H ) are the same: 
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In massive case we observe five independent HAs, while in the case of zero loop fermion mass (to / = 0) 
one gets only four independent HAs which are very compact: 
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All the others fcrmionic HAs arc equal to +1. 



There are relations among helicity HAs due to C,P,T-invariance. Moreover, there is another relation 
due to crossing symmetry: 



TL+ — + (s, t, u) = TL+-+- (s, u, t) , 
but this fact does not mean reducing the number of independent HAs. 



(15) 



3 Process 77 — > 77 in SANC 

3.1 SANC process tree 

In this section we briefly describe analytic modules relevant for 77—5-77. 

For boxes the SANC idea of precomputation becomes vitally important pQ . Calculation of some boxes 
for some particular processes takes so much time that an external user should refrain from repeating pre- 
computation in the SANC system. Furthermore, the richness of boxes requires a classification. Depending 
on the type of external lines (/ for fermions and b for bosons), we distinguish three large classes of boxes: 
////, ffbb and bbbb. 
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Figure 3: 77 — > 77 in the SANC process tree 
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The sum of contributions of fermionic loop boxes form a gauge-invariant and UV- finite subset, which 
is a consequence of Ward Identity Eq. (J9|. It is true for sum of bosonic contribution too. So we can 
distinguish QED and EW part of the process in the analitycal calculations. 

The precomputation file AAAA Box (see the SANC process tree in Figure [3] from SANC client [S]) contains 
the sequence of procedures for calculation of the covariant amplitude. At this step we suppose, that all 
momenta are incoming (denoted by p's) and photons are not on-mass-shell. Therefore, these results can 
be used for other processes which need these parts as building blocks. 

When we implement the process 77 — > 77 (see 4A QED or EW Processes branch), we use this building 
block several times by replacing incoming momenta p's by corresponding kincmatical momenta with the 
right signs, and calculate Ti by the module AA->AA (FF) , then helicity amplitudes by the module AA->AA 
(HA) and finally — the analytic expression for differential and total process cross section (XS) by the 
module AA->AA (XS) for the QED process. For the EW part we use numerical evaluation to get XS. 

3.2 Stand-alone SANC module 

Now let us introduce the SANC modules packages concept. 

In SANC system one has an opportunity of exporting the analytical results for numerical evolution [1 . 
Moreover, there are tools for checking the implementation of these Standard SANC FORTRAN modules 
(SSFM) j2] — the integrator of the process, based on the Vegas algorithm [15] . 

For light-by-light scattering process the SSFM are included in sanc_4A_vl.OO package, which is available 
for download from pj]. The numerical results are cross-section distributions from next section. Here we 
present the technical description of this package — some main flags and the options. 

I SANC Modules. 

I SANC NC package (Neutral Current processes modules). 

The previous versions of SANC NC package you can find in SANC Archives. 

O 29/07/2009 SANC NC vl.21 package (356 Kbtgz-file) [last stable version] 
In vl.21 package some bugs in the EW part are fixed. Details are in the file CHANGES. 

O 13/02/200B SANC NC vl.10 package (300 Kbtgz-file) [stable version] 

In vl.10 package the f lflHZ processes are added. For description of the flflHZ processes please refer to the 
paper D. Bardin et al., Comput. Phys. Comm. 177 (2007J 738-756. Details are in the file CHANGES. 

This package is intended for calculation of the 1-loop radiative correction to Drell-Yan Neutral Current processes at 
partonic level. Fora Technical Description of this module please refer to the paper Eur. Phys. J. C54 (2008) 451. 

I SANC CC package (Charged Current processes modules). 

The previous versions of SANC CC package you can find in SANC Archives. 

O 16/06/2009 SANC CC vl.30 package (236 Kbtgz-file) [last stable version] 
In vl.30 package the modules for new processes are added. Details are in the file CHANGES. 

© 05/12/200B SANC CCvl.20 package (160 Kbtgz-file) [stable version] 
In vl.20 package the modules for new processes are added. Details are in the file CHANGES. 

O 2B/10/200B SANC CC vl.ll package (131 Kbtgz-file) [stable version] 

In vl.ll package some bugs in the QCD soft-virtual part are fixed. Details are in the file CHANGES. 

This package is intended for calculation of the 1-loop radiative correction to Drell-Yan Charged Current processes 
at partonic level, see A. Arbuzov et al., Eur. Phys. J. C46 (2006) 407. 

I SANC 4A package. 

O 19/11/2009 SANC 4A vl.00 package (21 Kbtgz-file) [last stable version] 

This package is intended for calculation of the 1-loop full EW correction to light-by-light scattering processes, see 
D. Bardin et al.. "Light-by-light scattering in SANC". hep-ph/06111BB. 



Figure 4: SANC modules download web site [5] 
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The package files: 



• bbbb_main_ll_ll.F — the main file; 

• bbbb_ha_ll_ll.F — the HA from SANC system; 

• *_rl6.f — the library of special functions and algorithms; 

• *_input.h — the set of various setups of input parameters; 

• README, INSTALL and other instructions files. 

In README and INSTALL files one can find instructions how to use the package. The main options 
one can change in bbbb_main_ll_ll.F — the main file. 

There are some main flags and the options with next variants. 

• iha(I) — choice of helicities sum in cross-section: 
1=0 total helicities sum 

1=1 ++ helicities sum 
1=2 +— helicities sum 

• iqed(I) — choice of calculations for QED corrections: 
1=0 without QED corrections 

1=1 with QED corrections 

• iew(I) — choice of calculations for EW corrections: 
1=0 without EW corrections 

1=1 with EW corrections 

• gfscheme(I) — choice of the EW scheme: 
1=0 ao calculation scheme 

1=1 Gp scheme 

1=2 G'f scheme, a test option, when ao is replaced by aa F — \/2GfM^ (l — M^/M§) /it 

• isetup(I) — choice of the setup: 
1=0 Standard SANC input 

1=1 Les Houches Workshop (2005) 
1=2 Tcvatron-for-LHC Workshop (2006) 



1=3 Custom setup (17) 



• start (I) — choice of start point for 4 order logarithmic scale of yfs: 
I=lq-5 logarithmic scale y/s from 0.1 MeV to 1 GeV 
I=lq-1 logarithmic scale yfs from 1 GeV to 10 4 GeV 

To get full EW answer with interfarence one should set iqed=l and iew=l. 

3.3 Results and comparison 

To test our analytical results we calculate final answer for the total cross section for QED part in the 
massless limit after substitutions of helicity amplitudes and angular integration: 



' mt — 1 — + l 4- 2 ~ 8^ 2 C(3) + ^ 4 - 24 C (5) ) . (Ifi, 



77 _> 77 



2W 2 V 5 2 ,w 225 



This result was compared with [6] and the complete agreement was found. Also the limit of helicity 
amplitudes QED was compared separately with |llj and again full agreement was observed. The massive 
expression of HA for QED and E W parts were compared with [5] and [3] . 
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Figure 5: 77 — > 77 SANC and [H] cross-section (++ and H — ) comparison 
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The numerical evaluation of analytical expressions for QED, weak and full EW parts within SSFM 
sanc_4A_vl.OO package for cross section were comparied with jH] (see Figure [EJ. The setup is: cos 8 < 
cos7r/6 — cuts; My/ = 80.22 GeV; m t = 120 GeV; a = 1/137.036; other fermions are massless. 

Also according to [5] the sum of helicities cross-sections were calculated for MeV and GeV regions 
(see Figure [6]) with custom setup: 

• a =1/128; 

• cos0 < cos7r/6; 

• M w = 80.22 GeV; 

• m e = 0.51099892 MeV, 

• m u = 0.062 GeV, m c 

• m d = 0.083 GeV, m s 

4 Conclusions 

Let us summarise the results. 

In this paper the helicities amplitudes expressions for light-by-light scattering in general (massive) 
and in limiting (massless) cases were listed for QED (fermionic) and weak (bosonic) parts. 

We shortly described precomputation strategy of the SANC system [T] and the place of this process on 
the SANC process tree. 

The implementation of analitycal results and the SSFM concept were described. 

Additionally, we calculated and gave the answers for particular cases of D , C and B Q Passarino- 
Vcltman (PV) functions for massive and massless loop particles, strings and basis for covariant amplitude. 

The SSFM sanc_4A_vl.OO package was builded and tested. Its numerical results were comparied with 
those existing in the literature. The package is available for download at web pages 
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m M = 0.105658369 GeV, m T = 1.77699 GeV; 
= 1.50 GeV, rm = 120.0 GeV; 

= 0.215 GeV, m b = 4.7 GeV. (17) 
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Appendix 



To obtain the cross section of the 77 — > 77 process in an analytic form we have to compute in an 
explicit form the master integrals, Bq, Cq, Do — scalar PV functions [T2], [TJ — for a particular set of 
parameters. In Do and Co functions one can see collinear divergences, but the differential cross section 
is free of mass singularities which completely cancel in the sum of all terms. The Ao and Bq functions 
contain UV divergences, which cancel in the sum of box contributions. In the process of computation we 
face also a problem of the "angular edge" divergences, but they are not physical and cancel completely. 

I Bq function 

The Bq function for 77 — > 77 process reads: 



B (Q 2 ;M, M) 



In 



M2 



where 



and 



fl 2 



AM 2 



/?ln 



j3+l 
(3-1 



(18) 



M 2 = M 2 — ie . 



In zero limit of fermion mass: M — > 0, one gets: 

1 



B (s; M, M) 
B (t;M,M) 
B (u; M, M) 



l t + In 
lu + In 



(19) 



(20) 



II Cq function 

The Co function is: 

C (0,0,Q 2 ;M,M,M) 
After calculations: 



1 X 

J dx J dy (Q 2 y - Q 2 xy + M 2 ^j 







C (0,0,Q 2 ;M,M,M) = — 







Hi) 









(21) 



(22) 



where 



X lt2 = -(1±0). 



(23) 
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For M -> 0: 



C (0, 0, s; M, M, M) = 
C {0,0,t;M,M,M) = 
C o (0,0,u;M,M,M) = 



1 


In 


f M 2 






2s 




V s 


1 


'"( 








2t 






1 


In 








2u 







-i 2 



ITT 



(24) 



III Dq function 

The Dq function looks like: 

l x y 

D (0, 0, 0, 0, Q 2 ,P 2 ; M, M, M, M) = J dx J dy J dz ((Q 2 + P 2 ) xz - P 2 xy - Q 2 yz + P 2 y - P 2 z + M< 



After lengthy calculations: 



(25) 



where 



D Q (0, 0, 0, 0, Q 2 ,P 2 ; M, M, M,M)=- 



stA 3 y 



-Li 2 - 



A 



33 



Ai + A 3 



2 + 



- Li 2 [ - 



12 

A 3 3 



A 2 + A, 



-Lin - 



Li 2 

A 13 
1 + A3 



A33 
A 23 



- Li 2 - 



A 



23 



1 + A E 



^ ln2 (irk) 2C(2) + 2z7Tln i 1 + i) ^- /m ( Ai + - ( 26 ) 



/ 4M 2 
1 - , A 3 



AM 2 (Q 2 + P 2 ) 



'1- 



4M 2 (Q 2 + P 2 ) 
33 Q 2 P 2 (1 + A 3 ) ' 23 



4M 2 



Q 2 (A 2 + A 3 ) ' 



13 



Q 2 P 2 

AM 2 
P 2 {A l + A 3 y 



(x) is a function of Heaviside. 
For M -> 0: 

D (0, 0, 0, 0, s, t; M, M, M, M) ■■ 
Do (0, 0, 0, 0, u, t; M, M, M, M) ■■ 
Do (0, 0, 0, 0, s, it; M, M, M, M) = 



2 
si 
2 

ui 
2 



In 2 



hr 



In 2 



M 2 \ , / Af 2 \ 7T 2 . , 

M 2 \ , /M 2 \ , . tt 2 , , ■ 
— J - ln I — J (l t + l u ) - — + l t l u 



M 2 



+ ln 



M 2 
u 



IT 



+ in ln 



M 2 



M 2 
it 



(27) 



(28) 
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IV Strings and basis 



To obtain a compact form of structures of the amplitude we choose 14 auxiliary tensorial strings: 



a/3_ 
T l 


PlBP2a + 


\s8 a (3 , 


T 2 = 


PZ^PAv + 


±s5 


T 3 — 


P2uP38 + \t5f) v , 


• i 

T 4 = 


PlliPia + 






PlvPZa + 


2 LLU a.u •> 


T 6 — 


P4/3P2fi + \u5fin , 
Plfl - st- 1 p 3/3 , 


T 7 — 


Pi» - tu 


~ 1 P2 l i , 


T% = 


p\ v - ut~ 


1 P2u , 




T io — 


P2a - SU 


-1 P3a , 




Pifj, ) T 12 


= P3v, 


rfs = 


P2f}, T? A =p la . 



(29) 



The complete basis T"' ^ can be presented in a compact form with an aid of the auxiliary strings tj: 







afi uv 
T l '2 ! 


1 2 ~ 


Bv an 
T 3 T 4 ) 


rpapfiv 
X 3 ~ 


T 5 T 6 ' 


rpaBiiv 




ap n v 
T l T 7 T 8 i 


X 5 




V" P a 
~2 T 9 T 10> 


1 6 — 


u, a 
T 3 T 7 T 10> 


rpaPpv 
1 7 — 


an v P 
T 4 T 8 T 9 ' 


rpaBfiv 
1 8 




av H P 
~5 T 7 T 9 ' 


1 9 




^.Pp^v^a 
T 6 r 8 T 10> 


1 w — 


M v P a 
T 7 T 8 T 9 T 10' 


rpaBfi^ 

1 11 — 


aP u v 
T l T 7 T \2i 


rpaBfiv 
1 12 




^ P a 
T 2 T 9 T 14) 


1 13 




JPvjH„a 
T 3 T 7 T 14; 


rpO-PilV 
1 li — 


an v B 
T 4 T 8 r 13> 


rpaPpv 
J 15 " 


T 5 T 7'13' 


rpapn u 

1 w 




-P^-v-a 
T 6 r 8 T 14i 


1 17 




*P M i/ 
T l T 11 T 8 ' 


L \8 — 


T 2 T 13 T lCb 


L \9 " 


T 3 T ll'10> 


rpaBiiv 
1 29 




T 4 T 12 T 9 ' 


1 2l 




T ai/ T ^ 

T 5 T ll r 9 i 


rpOLp^lV 

1 22 — 


T 6 T 12 T 10' 


rpaBnv 
J 23 — 


T a P Li v 


rpaBfiv 
1 24 




-P-V-P -a 
~2 t 13 t 14j 


1 25 




T 3 T 11'14) 


rpaPitv 
1 26 — 


ap, v B 
~4 T 12'13' 


rpaBnv 
1 27 — 


T 5 T ll r 13' 


rpaBnv 
1 28 




T 6 r 12 r 14) 


rpafipv 

1 29 




~7 T 8 T 13 T 14 7 


rpaPixv 
J 30 — 


T 7 T 9 T 12 T 14; 


rpaBfj-v 
J 31 — 


v T P 
T 7 T 10 T 12 T 13' 


rpOpLlU 

1 32 




T 8 r 9 r ll r 14; 






i/ a /3 M 
T 8 T 10 T 13'11' 


rpOtfiflL* 

1 34 — 


P a fi v 
T 9 r 10 T ll T 12; 


rpaBn v 
J 35 — 


P 'a 
T 7 T 8 T 9 '14; 


rpapfi^ 
1 36 




A» ^ a /3 
T 7 '8 r 10'13' 


1 37 




M a _i/ 
T 7 T 9 r 10 T 12) 


1 38 ~ 


~i> P a M 

T 8 T 9 T io T in 


rpaPpv 
J 39 " 


M !v /3 a 
T ll r 12 r 13 T 10' 


rpaBiiv 
J 40 




_/* _y x a _/3 
T ll r 12 r 14 T 9 


1 4l 




_/* _/3 T a v 
T ll r 13 T 14 T 8 ; 


J 42 — 


T 12 T 13 T 14 T 7 ' 


rpaBnv 
J 43 — 


_M ^ _a 
T ll r 12 T 13 T 14- 









(30) 
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